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We have measured the dielectric constant and the pyroelectric current of orthorhombic (space group Pbam) 
NdCrTiOs polycrystalline samples. The dielectric constant and the pyroelectric current show features associated 
with ferroelectric transitions at the antiferromagnetic transition temperature (Tn = 21 K). The effect of magnetic 
fields is to enhance the features almost linearly up to the maximum measured field (7 T) with a spontaneous 
polarization value of ~ 3.5 /iC/m 2 . Two possible scenarios, the linear magnetoelectric effect and multiferroicity 
(antiferromagnetism + ferroelectricity), are discussed as possible explanations for the observations. 
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PACS numbers: 75.85.+1, 77.70.+a, 77.80.-e, 75.80.+q 



I. INTRODUCTION 

Multiferroicity, i.e. when two or more ferroic orders co- 
exist in a single phase of a material, has attracted renewed 
great interest during the last a decade. Especially, multifer- 
roics with magnetic and ferroelectric orders have been under 
intense investigation, both experimentally and theoretically, 
due to its possible application for memory devices. The com- 
mon features of this type of multiferroics include dielectric 
anomalies and the emergence of spontaneous electric polar- 
ization at the magnetic ordering temperature, which are also 
strongly dependent on the strength and orientation of the mag- 
netic field. In these materials, ferroelectricity is usually en- 
abled via different types of exchange interactions which break 
the spatial inversion symmetry of a paraelectric phase. On the 
other hand, magnetoelectric (ME) materials show strong mag- 
netic field dependent polarization and dielectric anomalies^ 
The consideration of magnetic symmetry and the thermody- 
namic potential can give a great deal of information to under- 
stand the ME effect. For example, both broken spatial- and 
time-inversion symmetries are required for a linear ME effect 
due to the linear potential term {EH) in the thermodynamic 
potential, while the spatial-inversion symmetry can be con- 
served for a quadratic term (EEH). 

Two mechanisms have been proposed as the macroscopic 
driving force for the symmetry breaking in a multifer- 
roic phase transition, namely, the spin-current models^ and 
magnetostriction 8 . Within the spin current model the sym- 
metry breaking is due to the Dzyaloshinskii-Moriya (DM) 
interaction^^ where the antisymmetric exchange allows an 
electric polarization, P, proportional to ey x Si x Sj, where 
eij is the unit vector connecting two sites of spin Si and Sj . 
The spin current is more relevant to a system with spiral mag- 
netic structure such as TbMnOyi, DyMnOy^, Ni 3 V 8 Og±2, 
CuFe02 14 and other compounds. 15 There are fewer examples 
to show multiferroicity with collinear magnetic structures. In 
Dy(Gd)Fe03^ and Ca 3 Co 2 - x Mi\ x Og^., the DM term van- 
ishes from symmetry and the magnetostriction mechanism 
was proposed to explain the multiferroicity. Another Pbam 
multiferroic irM^Os (R = Y, Ho, Tb, Er) has more compli- 
cated spin structures with orderings of different commensura- 
bility. Its origin of ferroelectricity has been ascribed to either 



the spin current mode l * or magnetostrictionisiiL Concern- 
ing the ME effect, regardless of the apparently different under- 
lying mechanism, these types of multiferroics generally show 
large ME effect, which sometimes is linear Pi = otijHj or 
Mi = OLijEj depending on the magnetic and crystallographic 
symmetries. 

The crystalline and magnetic structure of powder samples 
of NdCrTiC<5 has been investigated via X-rays and neutron 
diffraction in 1970. 22 The spins of the Cr and Nd ions were 
found to be antiferromagnetically (AFM) ordered below 13 
K with all magnetic ions participating in the magnetic or- 
der. The spins of the Cr ions are collinearly oriented along 
the c-axis, while the Nd moments are ordered in the afc-plane. 
Subsequently, the ME effect in NdCrTiOs was discovered^ 
NdCrTiOs was one of the first ME materials possessing two 
distinct magnetic sublattices, namely Cr 3+ and Nd 3+ , known 
at that time. The ME susceptibility was measured 2 ^ on a pow- 
der pellet sample showing a Tn of 20.5 K, rather than 13 
Ki22 The magnetic point symmetry (mmm') allows a linear 
ME effect^ Hence, under an applied electric field a magne- 
tization is induced, Mi = UijEj, where Ojj showed a rapid 
increase at Tn = 20.5 K and a drop below 8 The magni- 
tude of a.ij observed was of the order of 1 x 10~ 5 CGS units, 
i.e. smaller by one order of magnitude when compared to the 
well known ME material Q-2O3; 25 In Ref. |23|it is argued that 
below Tn the magnetic moments of the Cr 3+ ions order coop- 
eratively and then induce the order of the Nd 3+ ions by means 
of an exchange coupling, in contrast to a direct cooperative or- 
dering of the Nd 3+ spin subsystem. 

Assuming that the AFM order in NdCrTiOs is driven by the 
Cr 3+ spins, their collinear order would suggest that the mech- 
anism giving rise to the ME is magnetostriction. However, 
if both magnetic sublattices participate in the built-up of the 
cooperative AFM order, then the order has to be considered 
noncollinear and the spin-current model would be applicable. 
The open question of the origin of the coexistence of AFM 
and the ME effect motivated us to present a detailed study of 
the ME behavior of this compound. 

In this work, we measure the electric response (dielectric 
constant, polarization) of NdCrTiOs under magnetic field, 
which leads to observations of signatures of multiferroic and 
ME behavior such as a dielectric anomaly and spontaneous 
switchable polarization below Tjy. For the measured mag- 



netic field range up to 7 T, the saturated polarization (P sa t) 
increases with magnetic fields with a maximum P sat ~ 
3.5 /iC/m 2 . The observed ME annealing effect tends to point 
NdCrTiOs as ME material rather than multiferroic, but we 
cannot rule out the possibility that it is a multiferroic. 



II. EXPERIMENTAL 

Polycrystalline NdCrTiOs was made by solid state reac- 
tion. A stoichiometric mixture of Nd2C<3, Q2O3, and Ti02 
was ground together and calcinated in air at 1000 °C for 24 
hours. The powder was then reground and pressed into pel- 
lets under 400 atm hydrostatic pressure, and again calcinated 
in air at 1300 °C for another 24 hours. The x-ray powder 
diffraction (XRD) patterns were recorded by a HUBER Imag- 
ing Plate Guinier Camera 670 with Cu K a \ radiation ( 1 .54059 
A) with a Ge monochromator. The DC magnetic susceptibility 
measurements were made with a Quantum Design supercon- 
ducting interference device (SQUID) magnetometer with ap- 
plied magnetic field of H = 0.1 T. The specific heat measure- 
ments were performed in a PPMS (Physical Property Mea- 
surement System by Quantum Design). A sample of cylin- 
drical shape was cut from the pressed pellet with dimensions 
of 16.0 mm 2 x 1.0 mm for the measurements of the electrical 
properties and electrodes were painted with silver paste on the 
plane surfaces. The dielectric constant was measured using 
an Andeen-Hagerling AH-2700A capacitance bridge operat- 
ing at a frequency of 10 kHz. The pyroelectric current was 
measured using a Keithley 6517A electrometer on warming 
after poling the crystal in an electric field while cooling down 
from above Tn. The sample was short-circuited for 15 min- 
utes before warming up to remove any residual charges. The 
spontaneous polarization was obtained by integration of the 
pyroelectric current with respect to time. 



III. RESULTS 

The room-temperature X-ray diffraction (XRD) pattern 
of NdCrTiOs shows a single phase with an orthorhombic, 
Pbam, structure. The lattice parameters are a = 7.5812(2) 
A, b = 8.6803(2) A and c = 5.8123(3) A, similar to the pre- 
viously reported data, 22 The X-ray and neutron diffraction 
studvi 2 . identified that 95% of Cr 3+ and 5% Ti 4+ are dis- 
tributed over the 4f sites (center of oxygen coordinated octa- 
hedral) and 5%Cr 3+ - 95% Ti 4+ over the 4h sites (base center 
of a square pyramid) while Nd 3+ ions are located at the 4g 
sites. 

The DC magnetic susceptibility (see Fig. 2(a)) shows a 
slope change (kink) at 21 K. At the same temperature, the 
specific heat shows a A shape peak, which indicates that this 
transition involves long range magnetic order. This transi- 
tion temperature is in good agreement with that previously re- 
ported from ME measurements The earlier reported neutron 
diffraction study 22 has identified the low temperature phase 
(T < 13 K) as having AFM order. The magnetic moments 
of Cr 3+ in the 4f sites (Cr (4f)) are pointing along the c-axis, 
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FIG. 1: (Color online) The crystallographic and magnetic structure 
of NdCrTiOs according to Ref . |22|. Here we ignored the small inter- 
change of sites found for Cr and Ti atoms (see text). 



with alternating up-ward and down-ward orientations. The 
magnetic moments of Nd 3+ and Cr 3+ in the 4h sites (Cr (4h)) 
lie in the ab plane. The Cr (4h) moments are aligned along 
the &-axis and AFM correlated in that direction, while they 
are ferromagnetically correlated in the a-direction. The Nd 
moments, on the other hand, order in a similar way except 
that they are tilted away from the 6-axis by 12 dcg. The long 
range order of the Nd moments is believed to follow the or- 
dered Cr moments rather than being induced by a coopera- 
tive interaction of Nd 3+ spin systemsi 2 ^ 6 . The spin structure 
at temperatures well below Tn (ignoring the 5% of switched 
Ti and Cr ions between the 4h and 4f sites) and below the 
temperature of the magnetic susceptibility peak (T = 1 1 K) is 
depicted in Fig. 1. Note that the slope of the susceptibility 
begins to change below 13 K, which is reported in Ref. |22| as 
the temperature below which the magnetic Bragg peaks were 
observed. 

Fig. 3 shows the dielectric constant measured for different 
applied magnetic fields. For H = T, the real part of the 
dielectric constant (e') shows a slope change around Tn = 21 
K. When a magnetic field is applied, a peak is induced in e' as 
shown in Fig. 3(a). Since Tn decreases with magnetic field in 
an antiferromagnet (the order is gradually suppressed), with 
increasing magnetic field, the position of the peak moves to 
lower temperatures. Since the spin symmetry is broken by 
the applied magnetic field, the peak becomes stronger with 
increasing field. We also notice that (i) the position of the peak 
in e' shows no frequency dependence (not shown here); (ii) 
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FIG. 2: Temperature dependence of (a) the magnetic susceptibility 
and (b) the specific heat of NdCrTiOs. Each inset shows data over 
an extended temperature range around Tn . Note the kink in the sus- 
ceptibility at 21 K. 



the imaginary part of the dielectric constant, or the dissipation, 
also manifests a peak around a similar temperature as the peak 
of e' with applied magnetic field. These two behaviors prove 
that the anomaly of e 1 is really related to the phase transition. 
Hence, the observation of the anomaly of e' around the AFM 
ordering transition at Tn = 21 K demonstrates the coupling 
between the magnetic and electric properties of NdCrTiO^. 
We did not observe any anomaly associated with the broad 
peak in the magnetic susceptibility around T= 1 1 K. 

Fig. 4(a) shows the pyroelectric current measured for dif- 
ferent magnetic fields and a fixed electric field. The magnetic 
fields were applied during both cooling and warming up. In 
zero magnetic field, no pyroelectric current is observed. When 
a magnetic field is applied, a pyroelectric current arises in 
a temperature interval related to the phase transition. This 
broad peak of pyroelectric current starts at about 18 K and 
ends around 22 K and becomes stronger with increasing ap- 
plied magnetic field. Note that it is necessary to cool the sam- 
ple through Tn with an applied electric field to observe a py- 
roelectric current. The spontaneous polarization is obtained 
by integration of the pyroelectric current with respect to time 
and shown in Fig. 4(b). With increasing applied magnetic 
field, the polarization monotonically increases and saturates 
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FIG. 3: (Color online) The temperature dependence of the real part 
(a) and the imaginary part (b) of the dielectric constant of NdCrTiOs 
at different applied magnetic fields. The upper inset of (a) shows the 
dielectric constant with zero magnetic field over an extended tem- 
perature range around Tn. Note the kink at the antiferromagnetic 
transition. The lower inset of (a) shows the dielectric constant down 
to 6.5 K at 7 T. 



around 3.2 /iC/m for H around 6 ~ 7 T. The pyroelectric cur- 
rent and the onset of polarization are then directly related to 
the magnetic phase transition. Note that again no change in 
the pyroelectric current was observed around T= 1 1 K, where 
the magnetic susceptibility has its maximum. 

Fig. 5 displays the pyroelectric current and its integrated 
polarization measured on a sample from a different batch for 
various applied electric fields under a fixed magnetic field H 
= 6.7 T. The peak of the pyroelectric current around Tn and its 
related polarization become stronger with increasing applied 
electric field. Moreover, the sign of the pyroelectric current 
and polarization switches simultaneously with the sign of the 
applied electric field. As shown by the solid and open circles 
in Fig. 5, a similar pyroelectric current of opposite sign is ob- 
tained and hence a polarization of the same magnitude. This is 
typical for ME annealing conditions, i.e., poling with both H- 
and T-fields. The polarization increases linearly with applied 
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FIG. 4: (Color online) Temperature dependence of the pyroelectric 
current (a) and polarization (b) of NdCrTiOs for different applied 
magnetic fields for a fixed electric field of 580 kV/m (E_LH). The 
inset of (b) shows the magnetic field dependence of the low T polar- 
ization. 



FIG. 5: (Color online) Temperature dependence of the pyroelectric 
current (a) and polarization (b) of NdCrTiOs at different electric 
fields and a magnetic field H of 6.7 T (E_LH). The solid and open 
circles correspond to reversed electric fields. The inset shows the 
polarization as a function of poling electric fields. 



electric fields up to 200 kV/m followed by saturation behavior 
at higher electric fields. 

We studied the ME annealing effect more closely on an- 
other sample from the same batch measured in Fig. 5. For this 
study, we empolyed five different condition of ME annealings, 
namely, four different conditions of positive or negative com- 
bination of magnetic or electric fields and zero magnetic field 
cooling (ZMFC). As stated above, to observe polarization, the 
electric field is always applied during the annealing (cooling) 
process, and so is the magnetic field during the pyrolelectric 
current measurement (warming up). Fig. 6 shows the pyro- 
electric current data from this study when the electric field is 
applied perpendicular to the magnetic field. From this study, 
we found that the usual ME annealing condition (both mag- 
netic and electric fields applied during cooling) is necessary 
to induce polarization. The polarization reversal is achieved 
by changing the polarity of the electric field but not the mag- 
netic field. The same measurements were repeated after the 
sample was rotated by 90 degree so that the electric and the 
magnetic fields are parallel each other. The two different con- 
figurations (E_LH or E |jH) gave same results but with 50% 
larger polarization value for the E±H case. 



IV. DISCUSSION 

The maximum spontaneous polarization value (3.5 fiC/m) 
is smaller by two orders of magnitude as compared to other 
well-known multiferroic materials with the same structure, 
e.g., RMn2C>5 (R: rare earth elements or Y). The polarization 
value is expected to be larger in NdCrTiOs when measured 
for a single crystal, rather than for a polycrystalline sample, 
but is likely to remain within the same order of magnitude. 
Considering the almost linear increase of the polarization with 
magnetic fields (see the inset of Fig. 4(b)), larger magnetic 
fields could also enhance the polarization. From the linear re- 
lation, we obtained ME coefficients a xy = 1.5 x 10~ 5 CGS 
unit and a zz = 0.7 x 10~ 5 CGS units. These values are in the 
same order with what was obtained from the ME susceptibil- 
ity measurement^ 

Our findings, the dielectric constant anomaly and the po- 
larization reversal, suggest two possible scenarios; (1) multi- 
ferroicity with a concurrent AFM/FE transition or (2) a linear 
ME effect with AFM fluctuation and domain effects. A high 
resolution structural measurement or a polarization-electric 
field (PE) hysteresis experiment at low temperatures could 
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FIG. 6: (Color online) Temperature dependence of the pyroelectric 
current NdCrTiOs under different ME annealing conditions when the 
electric field is perpendicular to the magnetic field (E_LH). 

discern the two scenarios, but neither of them is available to 
us at the moment. At room temperature, NdCrTiOs is or- 
thorhombic (Pbam), a centrosymmetric structure at room tem- 
perature, which does not allow ferroelectricity via the DM in- 
teraction. Nevertheless, we cannot exclude the AFM/FE sce- 
nario without further definitive experiments. Hereafter, we 
discuss each scenario and its implications. 

A. Multiferroicity with concurrent AFM/FE transition 

In this case, the mechanism leading to FE is associated with 
the AFM phase transition. Due to the rather small magnitude 
of the polarization, the corresponding driving force inducing 
the symmetry breaking must be weak. To find the underlying 
mechanism, we examine the applicability of the spin current 
model first. The origin of the downturn of the magnetic sus- 
ceptibility below 1 1 K is not yet established. A possible plau- 
sible scenario is that the Nd 3+ spins (the J = 9/2 multiplet 
is split by crystalline fields to yield a Kramers ground dou- 
blet) gradually start to participate in the AFM order below 1 3 
K, a temperature reported in Ref. |22| as the onset of some of 
the magnetic Bragg reflections. The onset of the order of the 
Cr 3+ moments could not be determined in Ref. |22| and is seen 
as a very weak cusp in the magnetic susceptibility (see inset 
in Fig. 2(a)). At Xjv the fluctuations of the magnetic moments 
are rather weak compared to the thermal fluctuations at the 
transition. This picture agrees with the arguments presented 
in the earlier papers about this compoundi22i2^ Assuming that 
close to Tjv only the Cr 3+ moments order, then the magnetic 
structure is collinear, so that Si x Sj terms are zero in the tem- 



perature interval around Tjy. The spin current mechanism 6 7 
can then be excluded, leaving magnetostriction 8 as the candi- 
date to induce the FE. 

The above arguments are based on the assumption that the 
Nd 3+ and the Cr 3+ (4h) moments are not participating in the 
magnetic order until lower temperatures. Should this not be 
the case, then spin cross-products between Cr (4f) and Nd 
sites are nonzero and since the system lacks inversion sym- 
metry, the DM mechanism is active to produce a spin current. 
The strength of the spin current would depend on the magni- 
tude of the magnetic order in the Nd sublattice. Such mecha- 
nism, however, seems to be absent in NdCrTiOs considering 
that no features of the dielectric constant and the polarization 
were observed around the Nd ordering temperature. 

In order to explain the dielectric anomaly at zero field (see 
the inset of Fig. 2(a)) and the almost linear magnetic field de- 
pendence of the polarization we turn to the Ginzburg-Landau 
expansion of the free energy in terms of the electric and mag- 
netic fields.— Differentiation with respect to the fields leads to 
the polarization and the magnetization, 

P, (E, H) = Pf + eodj Ej + aijHj + ■ ■ • , ( 1 ) 
Mi{E,H) = Mf + f x 0lMj H j +a ij E j + --- , (2) 

where P s and M s denote the spontaneous polarization and 
magnetization, whereas e and fi are the electric and magnetic 
susceptibilities. The tensor a represents the induction of po- 
larization by a magnetic field or of magnetization by an elec- 
tric field which is known as the linear ME effect. For an AFM 
we have P s = M s = 0, so that in zero magnetic field an 
electric field can induce a finite magnetization, M, = ctijEj, 
and vice versa, in zero electric field a magnetic field can in- 
duce a finite polarization, Pj = otijHj. We have mentioned 
this effect in the Introduction. 

In Fig. 4(b) we presented the polarization as a function of 
temperature for several fields. The magnetostriction, when the 
sample is driven through the Neel transition, induces strains 
and hence a local magnetization in the grains and their bound- 
aries. The small displacements of the ions from their equi- 
librium position, the consequence of the magnetostriction, 
should change the polarizability of the sample and hence the 
dielectric constant. Also the local deviations from the ordered 
magnetization, #m,, arising from the magnetostriction should 
generate local electric fields through the a-tensor, giving rise 
to a change in the dielectric constant at TV. This effect is 
expected to be very weak in the absence of a magnetic field, 
but it is very magnetic field sensitive since all the local con- 
tributions Srrii align and become cooperative. This effect is 
restricted to the region around TV- It is also expected that the 
zero magnetic field anomaly of e' is stronger in a polycrys- 
talline sample, due to the numerous grain boundaries, than in a 
single crystal. Its magnetic field dependence should, however, 
be stronger in a single crystal. The effect on the absorption, 
i.e. e", is expected to be much weaker than for e' as indeed 
observed. 

The pyroelectric current generated when the system is 
cooled through the Neel transition is also expected to be pro- 
portional to SrriiEj. Hence, it is restricted to the region around 
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Tjv, grows with applied magnetic field approximately linearly 
and is linear in the electric field. The saturation polarization 
is then almost linear in the applied H, in agreement with the 
observations. The above is valid as long as the expression Pi 
= otijHj is valid, i.e. the fields are not too large; otherwise 
higher order terms in the Ginzburg-Landau expansion would 
have to be taken into account. 

A specific mechanism of how the magnetostriction ex- 
change can lead to a broken inversion symmetry in NdQTiOs 
below the Neel transition is not available yet. In other 
multiferroic materials with collinear spin structures, the in- 
teractions between neighboring spins from different sites 
seem to play a crucial role. For example, there are ferro- 
magnetically aligned Dy 3+ (Gd 3+ ) and Fe 3+ spin sheets in 
Dy(Gd)Fe03^ and ttU order of Co 2+ and Mn 4+ Ising 
spins in Ca3Co2- ;c Mn. ! ;06^ This kind of spin arrangement is 
not clearly apparent for NdCrTiOs according to the available 
spin structure (see Fig. 1). Similar measurements on a sin- 
gle crystal sample would give more insight on the preferred 
polarization direction and its magnetic field dependence. 

B. Non-FE transition/Linear ME effect 

First of all, it should be noted that the linear ME effect and 
the multiferroicity are not mutually exclusive, that is, they can 
occur in a single material as found in DyFe03. 16 Experimen- 
tally, it is easier to observe the linear ME effect compared to 
proving ferroelectricity. We attempted PE hysteresis mea- 
surement on NdCrTiOs samples but the experiment was in- 
conclusive since we could not observe typical hysteresis be- 
havior before the sample was damaged due to electrical break- 
down. 

There are characteristic behaviors observed in NdCrTiOs 
which favor the non-FE/linear ME effect scenario. For exam- 
ple, the ME annealing behavior shown in Fig. 6 is typical for 
powder ME material. It is well known that the ME annealing 
is necessary to observe ME effect in AFM material, otherwise 
random AFM domains result in zero net polarization!^— The 
reversal of the polarization was found in other ME materials 



when the ME annealing was done with opposite polarity of 
electric field^ 

Likewise, other experimental results can be also explained 
without evoking ferroelectricity; e.g. the linear field depen- 
dent polarization and the larger ME coefficient for (E_LH). 
The observed dielectric anomaly around TV is often found 
in other AFM systems, where spin fluctuations are ascribed to 
the anomaly^. The recovery of the dielectric constant well 
below the transition and the increased anomaly with the mag- 
netic field can as well be qualitatively explained by this model. 

However we can not completely exclude the possibility of 
multiferroicity. For example, it is also conceivable that the 
ME annealing is required for magnetic field induced FE tran- 
sition. The poling process is to acquire a single domain state 
while the sample is cooled down through a transition (for ex- 
ample, paraelectric to ferroelectric), therefore both magnetic 
and electric fields are necessary under the AFM/FE scenario 
as well. Between the two scenarios, however, we tend to in- 
cline to the second one, the linear ME effect rather than multi- 
ferroics. It is the simpler explanation and it is consistent with 
the small polarization. 

In summary, we have observed a field induced dielectric 
constant anomaly and electric polarization at the AFM state 
of NdCrTiOs. The polarization increases with magnetic fields 
and is reversible by changing the polarity of the poling elec- 
tric field. We discussed multiferroicity and linear ME effect 
as a possible mechanism to explain the experimental results. 
While it is difficult to draw a definitive conclusion, the lin- 
ear ME effect is more likely the cause considering the small 
polarization. 
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